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1 Introduction and summary 

Hilbert series of a supersymmetric gauge theory is a partition function that contains 
information about the spectrum of operators and the moduli space of vacuum solu- 
tions, and has been used to study many properties of a large class of theories. An 
introduction to and review of the subject is, e.g., the PhD thesis (l]. In this paper 
we focus on a narrowly defined calculational problem: computing explicit results for 
the Hilbert series of M = 1 supersymmetric QCD with U(N C ) or SU(N C ) color sym- 
metry and Nf flavors in the fundamental representation of U(Nj) flavor symmetry. 
The standard method to calculate the Hilbert series is based on the plethystic pro- 
gramme [2 13 and the Molien-Weyl formula 13 15 . This method typically leads 



to expressions with a large number of contour integrals, which makes it difficult to 
obtain explicit exact results or approximations for the Hilbert series. Recently, an 
alternative method was presented in fl6) (see also 17 for the SO(N c ) and Sp(N c ) 



symmetries), based on rewriting the Hilbert series as a basic object of interest in ran- 
dom matrix theory - a determinant of a Toeplitz matrix of Fourier coefficients of a 
periodic function, or in short a Toeplitz determinant. In random matrix theory they 
result from the problem of computing expectation values in an ensemble of unitary 
random matrices. With this method, a new set of tools, developed in recent years 
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in random matrix theory, become available and yield explicit results for the Hilbert 
series of SQCD. 

In this paper we discuss two additional techniques and calculate with them new 
explicit results for the Hilbert series. Our motivation comes from previous work where 
we have been applying and developing the techniques in a different context, which we 
now review briefly. Similar random matrix theory integrals or Toeplitz determinants 
as in 16 appear in the problem of computing amplitudes for the decay of unstable 
D-branes into various closed string or open string channels. In that context matrix 
integrals appear in an approach involving correlation functions in timelike boundary 



Liouville theory (TBL) [18], a scalar field theory with the action 

S = -— [ dzdz <90<90 + X f dr e^ (r) , (1) 

where the last term is a deformation by a Liouville term inserted at the boundary of 
the worldsheet, and the kinetic term of the scalar has a negative sign. A straight- 
forward path to compute correlation functions is to treat the boundary deformation 
term as a perturbation and expand it as a series in the path integral. After perform- 
ing the contractions, the iVth term in the perturbation series can be recognized as 
a U (N) matrix integral, with the standard Haar measure coming from the contrac- 
tions among the exponential terms [l9-2l].^] For explicit results for the correlation 



functions, one needs to compute the U(N) integrals. This becomes more difficult for 



n-point correlation functions with n > 2 [23] 24 , and involves similar integrals as in 
the Hilbert series problem: 



N , N n 



Z n ({w t };N) = ^y/n^ II \*** ~ II II I 1 " > ( 2 ) 

a=l l<a<b<N a=l i=l 

where Wi are n points in a unit disk in complex plane. We will review the Hilbert 
series and its connection to (J2j) in Section [2] 

We discuss two techniques that have been applied to the integrals Q, and will 
hence also apply to the Hilbert series: i) the use of orthogonal polynomials on the unit 
circle (Schur polynomials), and ii) relating the integrals to partition functions of log- 
gases with external charges, and then using electrostatics for saddle point estimates. 
Schur polynomials have appeared in the context of Hilbert series as the characteristics 



of SU{N) representations 13 16 and they have also been used extensively in the 
analysis of the M = 4 super Yang-Mills |25]|^] In this paper we derive several new 
results for the Hilbert series based on standard properties of Schur polynomials. The 
second technique to be discussed in this paper, the log-gas approach, has been studied 



by the authors previously in 24 ,27 as a way to obtain useful approximations to 
integrals 



1 Incidentally, as a precursor of connections to Yang-Mills theory, the appearance of U(N) matrix 
integrals inspired [22] to speculate that TBL is related to QCD in two dimensions. 

2 It is interesting to note that recent studies of Af = 2 dualities utilize Macdonald polynomials, 



which are generalizations of Schur polynomials 26 . 
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The paper is organized as follows. We begin by very briefly introducing the 
refined and unrefined Hilbert series, followed by a brief review of Schur polynomials. 
Then, using Schur polynomials, we derive an explicit formula for the refined Hilbert 
series when Nf < iV c . The formula is a summed simplification of a series expression 
conjectured in 16 . As a special case, the formula reduces to the result of 16 for 
the unrefined series. When Nf > N c , for the refined series we obtain an expression in 
terms of sums of Young tableaux. We also derive explicit results in a few cases with 
small Nf — N c . We then analyze the singularity structure of the (unrefined) Hilbert 
series. A direct calculation verifies that the orders of the poles of the Hilbert series 
agree with the dimensions of the moduli space of SQCD also for Nf > N c , where 
only a few explicit results are known. For a general discussion on the classical and 



quantum moduli spaces of vacua of supersymmetric gauge theories, see 28 



We then move to the second approach, the use of (electrostatics of) log-gas of point 
particles confined to the unit circle (also called the Dyson gas [29] ). We rewrite the 
unrefined Hilbert series as a canonical partition function of a log-gas (with the num- 
ber of particles depending on N c ) with an external charge (with the charge depending 
on Nf). This allows us to observe a physical phenomenon that is associated with the 
Nf = N c threshold. In the Veneziano limit of large Nf, N c with fixed ratio Nf/N c , it 
is possible to obtain an asymptotic approximation for the Hilbert series, by approx- 
imating the log-gas of charges with a continuous charge distribution. The leading 
contribution comes from the electrostatic minimum energy configuration for the con- 
tinuous charge distribution interacting with the external charge. When Nf < N c , the 
continuous distribution has nonzero density everywhere, and the leading approxima- 
tion for the Hilbert series agrees with the result from the Schur polynomial method. 
When Nf > N c , exact methods to compute the series lead to more complicated 
results. In the Veneziano limit, in the electrostatic method we observe that a gap cre- 
ated by the external charge may appear in the continuous charge distribution]^] This 
is reminiscent of 



32 



and we also find a third order phase transition. The asymp- 
totic result for the Hilbert series is still straightforward to find, and it again is simple 
to read off the dimension of the moduli space of SQCD. We also compare our ap- 



proach to that of 16 based on the Toeplitz determinants and evaluating them with 



the Geronimo-Case-Borodin-Okounkov (GCBO) formula. We show that this method 
works in the region where no gap appears in the asymptotic charge distribution of 
the log-gas method. 

Let us summarize in the following the main new results of this paper: 

• Color group U(N C ) 



We prove the general expression (15) for the refined U(N C ) Hilbert series 
as the sum over semi-simple Young tableaux in Section |2.3 



We obtain explicit formulas for the refined U(N C ) Hilbert series for Nf < 



N c , N f = N c +1, and N f = N c +2, in Eqs, (16), (19), and (20), respectively. 



3 A similar observation has been made in mathematics literature 30 . Another, earlier related 



result is 31 which studied the phase structure of unitary matrix models with polynomial potentials. 



A difference is that in 31 the electric fields which arise are not created by point-like external charges. 
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For the unrefined series, we work out the order of the (only) non-trivial 



pole in Eq. (34) which matches the dimension of the moduli space. 



We calculate the unrefined Hilbert series, and the residue of its pole, in 
the Veneziano limit up to corrections suppressed by 1/N% for all reasonable 
values of the parameters. Results are given in Eqs. (49), (56), and (57). 



Color group SU (N c 



We prove the general expression (24) for the refined SU(N C ) Hilbert series 
as the sum over semi-simple Young tableaux in Section |2.4 

We obtain explicit formulas for the refined SU(N C ) Hilbert series for Nf < 
N c , N f = N c , and N f = N c + 1, in Eqs, p6b, (pb, and (|30b, respectively. 



We work out the orders of most of the poles of the unrefined series in 



Eqs. (38) and (39), and find again agreement with the dimension of the 



moduli space. 

We prove that the results for the unrefined U(N C ) Hilbert series in the 
Veneziano limit also apply to the SU(N C ) Hilbert series, with certain re- 



strictions on the values of its parameters (see Eq. (61)). 



2 Hilbert series and Schur polynomials 



2.1 Hilbert series 



As in 16 , we begin by considering an M = 1 supersymmetric U (N c ) Yang-Mills the- 
ory with Nf flavors of quarks and anti-quarks in the fundamental and antifundamental 
representation. The theory has various global symmetries; of particular interest here 
are the U(1)q (U(1)q) charge symmetries, with quarks (anti-quarks) carrying charge 
+1 (—1), and the SU(Nf)i x SU(Nf) 2 flavor symmetry. The refined Hilbert series 
is a partition function jTJ,[Toj involving fugacities associated with the color and fla- 
vor indices and with the global charges. We denote the SU(Nf)i x SU(Nf)2 flavor 
fugacities by Xi, yi, with i = 1, . . . , Nf and the U(N C ) color fugacities by z a = e lTa , 
with a — 1, ... , N c . The U(1)q x U{1)q fugacities are denoted by t, t. With the help 
of the Molien-Weyl formula from representation theory, the refined Hilbert series can 
be written in a simple-looking form which however involves nested integrals. After 
rescaling the flavor fugacities, introducing 



x 



t- ~ - , i x 2 x 3 1 

I Xi, #2, • • • , XN f ) — \X\, , , • • • , 

X Xi x 2 X Nf -l 



~ ( ~ ~ ~ \ /I 2/12/2 
y = (2/i,J/2, • • • ,VN f ) = - , 3/jvv-i 

V2/i 2/2 2/3 

the refined Hilbert series becomes a Selberg type integral 



(3) 
(4) 



*w.) (*■ *> **v) = w ft P ^ I A W I 2 nfln 

c ' 0=1 J0 1=1 0=1 ^ 



tf = l (1 ~tXiZ a X ) (1-tyiZa) 



(5) 
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where A(z) is the Vandermonde determinant 



II ( Z a~Z b ). (6) 
l<a<ft<iV c 

Notice that 

N f N f 
j=l i=l 

For the reduction to SU(N C ) color symmetry, we introduce the extra constraint 
that the phases r a sum to zero (modulo 2tt): 



c a =l' /u fc=-oo \ a / 



9N f ,su(N e )(t,t,x,y) = — jj / — \ ( 2^ r„ - 2m A- J |A(s) 

x nn7r 



i=l a=l 



tXiZ a X ) (l - %Z a ) 



As pointed out in 16 , this result can be written as a sum over unconstrained integrals 

as 



9N f ,su(N c )(t,t,x,y) = ^2 h(t,i,x,y) 

k=—oo 

4(M>,y) = tttII/ ^|A(z)| 2 ^4---4 c 

iVc - a=i ^ Z7r 

N c 

x nn (1 _ % -i )(1 _^ a) ■ ( 9 ) 



j=l a=l 



where J equals the U(N C ) partition function gN s ,u{N c )- 

The unrefined Hilbert series are obtained by the setting x\ — • • ■ = x^ f = y\ — 
• • • — VNf = 1 in the above definitions, and we shall denote them by gN f ,u(N c )(t,i) 
and by 9Nf,su(N c ){t,t) for the U(N C ) and SU(N C ) gauge groups, respectively. 

It is well known that the matrix integrals ^ and ^ can be rewritten as N c x N c 
Toeplitz determinants [16). For the case of unrefined series, the elements of the deter- 
minants can be evaluated explicitly (see Appendix [A]) , leading to explicit expressions 
for the unrefined series for small N c and arbitrary Nf. While Toeplitz determinants 
can also be defined for the refined series, it is more efficient to use the language of 
Schur polynomials, which we shall discuss next. 

2.2 Schur polynomials 

The integral ^ has been previously analyzed in the context of D-brane decay 22 



27 . The method presented in [22 1 uses expansions in a particular set of symmetric 
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polynomials called the Schur polynomials, and applies directly to the present case. 
For completeness, we shall review the calculation here. 

To begin with, let A = (Ai, A2, • • • , A n ), be a partition of |A| = J2i with Ai > 
A 2 > ■ ■ ■ > A n > 0. Then A parametrizes a Young tableau, corresponding to a 
Young diagram where the ith column contains Aj boxes. We also define the length 
of the tableau £(X) = max{i : Aj > 0}. The Schur polynomials of the variables 
z — (z 1 ,Z2,---, z n ) are characterized by A: 

det ( z ^-i + i+i-A det ( z K- i+ i+i-i 



/ \ / \ _ / i,j=l,...,n _ \ / i,j=l,...,n 

s\{z — s\{Zi, z 2 ,..., z n — -. i n — — — , 

det (■-.; % , „ A ^ 

(10) 

where the denominator is the Vandermonde determinant. It is not difficult to see that 
this expression indeed defines a polynomial which is symmetric under the interchange 
of any two variables. 

We shall need the extension for the definition of the Schur polynomials to the cases 
where the length of the partition exceeds the number of variables, A = (Ai, A2, • • • , A m ) 
with m > n. It is convenient to define s\(z) = whenever £(X) > n, i.e., the Young 



diagram is larger than the number of variables, whereas the above definition (10) ap- 
plies for £{X) < n. With this extension, we have s\(zi, z 2 , z n , 0) = s\(zi, z 2 , z n ). 
The polynomials are orthogonal with respect to the Haar measure of 



1 n r 2n dr- n 
— r TT / -w 1 \A{.z)\ 2 s x (z 1) ... ) z n )s K (z 1) ... ) z n ) = T\6x k ,K k = 

i=i Jl) k=i 



as can be verified by direct computation using the definition {lO|)0 

The most important tool for us is the Cauchy identity. Let z = (zi, . . . , z n ) and 
w = (w±, . . . , w m ). Then 



n m 

i=l j=l 1 3 \ 



Here the sum over A is restricted by the requirement that both s\{z) and s\{w) are 
nonzero, which leads to £(X) < min(n, m)n 

The following property of the Schur polynomials, which follows easily from the 



definition (10), will also be useful. For k > 0, we have 

z k 1 z k 2 ---z k n s K (z) = s K+k (z) (k>0) , (13) 

where n + k is defined by (« + k)i = Ki + k with i — 1, . . . , n. We shall also denote 
k + k = k + kin) to stress that k boxes were added only to the first n columns of the 
Young diagram whenever this is not clear from the context. 



4 The orthogonality property holds for I (A), £(k) < n. If the length of either of the Young diagrams 
exceeds n, the integrand vanishes. 

5 Notice that if, e.g., £(X) < n, we take Ai = for n > i > £(X) in the definition of the Schur 
polynomial s\(z). 
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2.3 Results for the refined U(N C ) series 

Let us then apply these results to the refined Hilbert series of Eq. (J5j). Applying the 
Cauchy identity (with (n, m) = (Nf,N c )) we have 

N f N c 

n n (1 _ ii^j (i - tfec) = g ^^w^.ao'.w • (") 

where £(A),£(k) < min(Af,A c ) and we used s\(tz 1 , . . . ,tz Nf ) = t^s\(zi, . . . , z^X 
Using the orthogonality relation of Eq. ( pTj ), we get 

9Nf,u{N B )(t,i,x,y) = ^2t lxl s x (x)$ Kl s K (y)5 x , K 

A.K 

£ (ttT'^My) , (15) 

\:£(\)<imii(Nf,N c ) 

for all Nf,N c . If Ay < N c , the constraint £(X) < min(Ay-, A c ) = Aj is automatically 
satisfied. In this case we can use the Cauchy identity again, and get 

N f N f 

9N f ,u{N B ){t, t,x, y) = J] J] , (A/<A C ) (16) 

i=l j=l 1 

which is independent of A c . 

Ref. 16 conjectured the formula (15) for the refined Hilbert series, in the form of 
an infinite series with coefficients given by SU(Nf) characters (eqn. (2.79) in 16 ). 
We have rewritten the formula in terms of the Schur polynomials, and proven it here. 
For Nf < N c , the infinite series expansion then becomes equal to the right hand side of 
the Cauchy identity (12), so that the expansion reduces to the form (16) above. This 
explicit result for the refined Hilbert series is manifestly a generalization of Eq. (2.60) 
in 16 for the unrefined series, which was obtained using the Toeplitz determinant 



method and the Geronimo-Case-Borodin-Okounkov formula. 

Ref. 16 also contains explicit results for the unrefined series for Nf = N c + i with 
i = 1,2,3. Below we will work out the refined series for the two former cases; our 
results are consistent with those in (16) upon refining. 

In our case, when Nf > N c we have an extra constraint £(X) < N c which prevents 
us from doing the sum in a straightforward manner. However, small Nf — A c can be 
worked out. Let us start with Nf — N c + 1. We can write the refined series as a sum 
of two contributions: 



9N c +i,u(N c )(t,t,x,y) = (tiy Ms x( x ) s x(y) 

x-.e(\)<N c 

= ^(tt)l A l SA (£) SA (y) - i^sxWsxiy) . (17) 

A X:£(X)=N c +l 

The first term can be readily summed, and with the second one we proceed as follows. 
Since £(X) — Nf — N c + 1, all column heights Aj > 1, with i — 1, . . . , Nf. Therefore 
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there exists a partition k such that k + 1 = A, which by definition means that 
Ki + 1 = Aj for all i — 1, . . . , iV/. For all partitions k, the shifted partition k + 1 is 
included in the sum of the second term since £(k + 1) = Nf. We also notice that 
s K+ i(x) = (Y[i%i) s n(x) = s K (x) by Eqs. (jTl) and (13). Collecting these observations, 
we get 

(tt)Ws,(x)s x (y) = ]T(^ +1 Wi(^ + i(y) 

X:£(X)=N c +l K 

N a +1 



i,j=i - — 



Therefore we find 

iV c +l J 

9n c +i,u{Nc) (*> 2/) = (1 " («T C+1 ) I] i _ „- - • ( 19 ) 

ij=l 1 UX iVj 

The case Nf = N c + 2 can also be treated in a similar manner, but some extra 
tricks are necessary. The calculation is done in Appendix [Bj We find 



H v +->.(( \' e )(t,i,x,y) = Yl 1 _ t i 5 ,.~. 



X 



iVi 

n£(i - nfii(i - ttstiiu) 



It may be possible to simplify the sum in the square brackets further. One can also 
derive similar results for higher fixed Nf — N c , but the expressions quickly become 
very cumbersome. 

2.4 Results for the refined SU(N C ) series 

The above analysis can be repeated to a large extent in the case of SU(N C ) SQCD. 
Now the integrand in Eq. ^ contains 

(21) 

We proceed by using the identity (13). For k > 0, we have 

z\z\ ■ ■ ■ z k Nc s R (z) = s K+k (z) (k > 0) , (22) 

where k + fc is defined by (k + k)i = Ki + k with i = 1, . . . , iV c . For k < we use 
instead 

zfz£ • • • z k Nc s x (z) = s x+lk{ (z) (k < 0) . (23) 
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By combining the above equations, and by using orthogonality we find 

oo 

9N f ,su(N c )(t,t,x,y) = h ( 24 ) 



k=— oo 



1 _ At- ' J " c 

±u ^iA W i^-4„nn(i 

c a=l J ° i=l a=l v 

A,/t 

^t |A| s A (5)t |K| s K (y)5 A+ | fe | (iVc)iK ; 



tXiZ' 1 ) (1 - %z a ) 



(25) 



k > , 
fc < . 



A,K 



Recall that the partitions must satisfy £(X),£(n) < mm(Nf, N c ). We can now analyze 
separately different cases: 



Nf < N c : In this case the Kronecker delta symbols in (24) can only be nonzero 
for k = Oj^] Consequently, the integrals I k vanish for k ^ 0, and the result for g 
is the same as for U(N C ) SQCD, 



N f N f 

9Nf,SU(N a )(t, 1 X, = J = J] II 1 _ fZz, ~ 

i=\ j=i 1 ux iyj 



(N f < N c ) . (26) 



t Nck I I I I — -~ — ; k > o 

i=l i=l 

^> A| s A (x)t' 

A 

t=i i=i 

where we recalled the definitions of x and y in terms of x and The sum over 
k can also be done explicitly, giving 

9N c ,su(N c ){t,t,x,y) 

1 



Nf = N c : Now the Kronecker delta symbols in (24) can be nonzero for all values 
of k. We get 

K K 

N c N c 

ni L w~ ; k -° (27) 

^U,^A[+^l*[ aA+|Jfc|( y) = FcW^I . . . y]*" ^tW aA (i)tW aA+w (y) 

A 

n n i _ t i~ x .y. 



fc=— oo 



+ 



-, iV c AT C 



i - (ti) 



N c N c Nc 



i=i j=i 



i 



(29) 



6 The partitions are constrained by £(X),£(k) < N f < N c , but £(k + k(N c )) = N c for k > and 
£(A + |fc|(JV c )) = N c for fe < 0. Therefore, the Kronecker deltas vanish, unless k = 0. 
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Setting here the x's and y's to unity we recover the unrefined series calculated 
in 



16 



Nf — N c + 1: This special case is worked out in Appendix [Bj We find 

N f 

9n c +i,su{n c )' ' 



(t, t, x,y)=l[ 1 - (30) 

^ (i - (tt) N i(x k y e )) nS(i - tfckVi) n£i(i - ttxM 



x ^ 



^^(i - t N °/x k )(i - t N -/y £ )(i - ttxkyt) n^fcOi - x k) Hj&Wj 



Nf > N c with generic Nf-. In this case all partitions in (24) have lengths less 
than equal to N c . We can do the sums over the partitions by using the Kronecker 
deltas: 

t kN c t^s K+HNc) (x)t^s K (y) ; k > , 

tWN. t\ x \ Sx (x)t\ x \ Sx+MNc) (y) ; k<0. 

\:£{\)<N C 

Proceeding further is difficult, and we will later find an interpretation by using 
the log-gas approach, with qualitatively different behavior separated by the 
N f = N c threshold. 

2.5 The asymptotics of the unrefined Hilbert series 

Let us then discuss the asymptotics of the unrefined Hilbert series and the singularities 
of their sums for U(N C ) and SU(N C ) theories. As the series for Nf < N c were solved 
explicitly above, their singularities are readily known, and we can restrict to Nf > N c . 



For the unrefined series, the results (|15|) , (24), and (31) may be written as 



9N f ,uiN c ){t,t) = (")' A|rf A (32) 

\:£(\)<N C 
oo 

9N f ,su(N c )(t,t) = s ^t kNc (ti)^d K+k{Nc) d K 

k=0 k:1{k)<N c 
-1 

fe=-oo \:£(\)<N C 

where d\ = s\(l, 1, . . . , 1) denotes the dimension of the SU(Nf) representation char- 
acterized by A. 

The asymptotics of the series and correspondingly the singularities of the series 
can be found by a "brute force" calculation, which is sketched in Appendix [C| We 
find that for large k and s — \X\, 

G(s, k)= E d x+k{Nc) d x ~ (s + tfQWt-N.)^,*.-! (33) 

A:]A]=s, £{\)<N C 



11 



Consequently, we get 

oo oo 

9N f ,u(N c) (t,t) = £>")-G(a,0) ~ ^(tt)'^--^- 1 ~ (1 _ t - )2JV/jVe . jVg (34) 



s=0 



for Afy > AT C as — > 1, where the order of the singularity matches with the expected 
dimension of the moduli space. 

In the SU (N c ) case, let us first analyze the two terms of Eq. (32) separately. That 
is, we first write 

oo oo —1 oo 

9N f ,su(N c) (t,t) = ^2^2t kNc (tt) s G(s, k) + ^^y ^ l fc D • ( 35 ) 

fc=0 s=0 k=—oo s=0 

Here the first term is singular at t — > 1 and at tt — > 1. In the first case we find 



oo oo 

. V,,,^ v V '- J 

fl _ t )N c (N f -N c )+l 
k=0 s=0 v ' 



i i 

t ~i (l _ t )JVc(JV/-Jv c )+i (1 _ ^jvjv, ( 36 ) 

where the first line gives the order of the pole at t — 1, and the second line adds the 
leading behavior of its residue as we further take tt — > 1. For the other singularity 
we find similarly 

oo oo 1 

yy t kN ^ti) s G(s,k) ~ - — 4— , 2 y^ 



fc=0 s=0 v ' fc 

1 1 



I 1 

(37) 



t-H (1 _ t t)2N c Nf-m I _ f 



so that the result depends on the order of limits. For the second term in Eq. (32) we 
find analogous results with t •<->■ t. 
Adding the results up, we find 

i 

9N f ,SU(N c )(t,t) 



t^l (1 _ t ^N c (N f -N c )+l 
9N f ,SU{N c )(t,t) _~ (1 _ ^ Nc ( Nf -N c )+l ( 38 ) 

where Nf > N c ^ We are also tempted to write down the result as tt — >■ 1. However, 
the corresponding singularity appears at the same order in both terms of Eq. (35), 
and since taking tt — > 1 fixes t = 1/t, cancellations between the two terms can, and 
typically do, take place. Therefore, we cannot predict the order of the singularity at 
tt — > 1, expect that it is smaller or equal to 2N c Nf — N^. 



7 Here we only wrote down the order of the pole and left out the structure of the residue, which 
is given in Eq. (36 1 as tt — > 1. 
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Finally, let us calculate the order of the "overall" singularity, taking first t = t 



and then t — > 1. In this case, the two terms in Eq. (35) are identical for \k\ > 1 and 
add up without the possibility of any cancellations. As seen from Eqs. (36) and (37) 
by setting t — t, or by direct calculation, 

9N f ,su(N c ){t,t = t) ~ x - _ ^J fNc _ N 2 +1 (39) 

so that the order of the pole is indeed the dimension of the SQCD moduli space. The 
extra factor of (1 — t) with respect to the U(N C ) result arises from the additional sum 
over k. 



3 Hilbert series in the Veneziano limit 

There are two motivations to consider the Veneziano limit N c , Nf ^> 1 with fixed 
ratio Nf/N c . On one hand we can expect to find simplifications for Nf > N c where 



only a few explicit results are known. This was also done in 16 and it is interesting 
to compare the results. On the other hand, we will use a different method where 
we reinterpret the Hilbert series as a canonical ensemble partition function of a two- 
dimensional gas of point particles confined on a unit circle with an additional external 
charge. The system is a specific case of a Coulomb gas; we will follow a random 
matrix theory literature convention and refer to it as a log-gas (due to the logarithmic 



Coulomb potential in two dimensions), see 33 34. This method will bring new 
physical insight to the threshold Nf = N c . 

3.1 The log-gas approach and the U(N C ) Hilbert series 

We begin by rewriting the U(N C ) unrefined Hilbert series as a log-gas partition func- 



tion. From the series expansion (15) we observe that gN f ,u(N c )(t,t,x,y) and thus also 
the unrefined series gN f ,u(N c )(t, t) depend on t and t only through the combination tt. 
Therefore, we can take t = t without loss of generality. The defining integral can be 
written as 

, N c „ 2 7T i N c i 



c a=l J0 a=l 



We restrict to the region < t < 1 where this integral and the Hilbert series converge. 
We can rewrite the integral in the form of a canonical partition function 



1 N c r 2n j 
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where the Hamiltonian H is 

h = - lo s \ eiTa - eiTb i + N f Yl log \^ Ta ~ * i ' ( 43 ) 

l<a<b<N c a=l 

and we fixed the inverse temperature to /3 = 2. The Hamiltonian describes an ensem- 
ble of N c point particles (of charge +1) on the circle at e lTl , . . . , e iTjVc , which interact 
via the Coulomb potential, with an external particle with charge —Nf, located on 



the real line, inside the unit disk at t E (0, 1) [20,35 . The system is symmetric with 
respect to an overall sign flip of the charges. 



We shall now apply the general method described in 24 , 36 ; the case with one 



external charge inside the disk was studied in detail in 27 . The same method was 
used to analyze unitary matrix models with different potentials in 31 . We repeat 
the main points. The leading contribution to log gN f ,u(N c )(t) in the Veneziano limit 
of large N c j with fixed Nf/N c = r may be calculated by going to the continuum 
limit of log-gas. Calculating the partition function in this limit resolves into solving a 
potential problem with the external charge at the point z = t. There is an interesting 
transition point [27] : because all particles on the circle carry the same charge, the 
external charge may cause a gap in the charge distribution on the circle if it comes 
sufficiently close to it. Therefore, we actually need to solve two potential problems: 
one with a gap and one without. In comparing the analysis with that in 27 , notice 
that because the sign of the charge of the external particle is now opposite to that 
of the distribution on the circle, their Coulomb interaction is attractive, and the gap 
will appear at the opposite side. 

As pointed out above, there will be a threshold for the gap formation at Nf = N c , 
which is qualitatively understood as follows. As the external charge approaches the 
boundary, it is able to attract at most Nf positive unit charges in its vicinity, because 
these charges will screen the external charge so that its potential effectively vanishes 
at long distances. In particular, if Nf > N c , all unit charges are attracted by the 
external charge and a gap eventually forms opposite to it. For N c > Nf, the leftover 
N c — Nf charges will remain evenly distributed on the unit circle, and no gap is 
formed. This behavior will be confirmed by the explicit calculation below. 

No gap in the distribution. The solution for a potential on the disk with one 
charge is well-known, 

V{z) = N f \og^ t . (44) 
The charge distribution on the boundary of the disk (with z = e tT ) is given, up to a 



constant, by the radial electric field 27 : 

p(x) = const + l\V'(,)\^ = 1(nc~N, + N f l + t l:* cixT ) . (45) 

Here the constant was fixed by the normalization § drp{r) = N c . The distribution is 
minimized opposite to the external particle, 

/ x / x N c ( 2tr \ 
minp(r)=p(7r) = ^h -— j , (46) 



14 



where r = Nf/N c . A gap forms when the minimum value reaches zero. The gapless 
phase requires 

K^-^U. (47) 
Note that for Nf < N c , \t c \ > 1 so that the system is always in the gapless phase, as 



expected. The continuum energy of the system is as in 27 



£ = ~\ J dndrz p(ri)p(r 2 ) log |e iri - e iT2 \+N f J dr p{r) log \e lT - t\ 

= -J-\ og (l-t 2 ) . (48) 

The leading result for the unrefined Hilbert series in the gapless phase reads 

9N f ,u { N c) {t)^e- 2£ = 1 ; < t < t c . (49) 

(1 -t 2 ) f 

The constraint for t can be dropped if Nf < N c , and the gapless result holds for all 
< t < 1. We immediately notice that the order of the pole at t — 1 agrees with 
the dimension of the moduli space of U (N c ) SQCD for Nf < N c , and actually the 
results matches with the exact result in this region (see Eq. (16) and fl6]). We shall 



also compare this to the Schur polynomial result and the results in [16] , but first we 
discuss what happens when a gap forms. 



One gap in the distribution. For t > t c the charge distribution develops a gap 
opposite to the external charge. For this to be possible, we need Nf > N c . The 
calculation of the electrostatic energy £ is quite technical in the presence of the gap. 
We only report the results here, as a similar calculation (for opposite bulk charge) 



was carried out in detail in 27 . See also 36 for a generic discussion of the solution 



for arbitrary numbers of gaps. 



Following 27 , we see that the support of the charge distribution is 



The 



edges of the gap, ±r c , are given by 



sin 



2 T c 



(r-1) 2 



(1+*) 2 
(l-t)a 



i m 2 - 1 ' 



(50) 



where 



X 



5(t) 



l + t 



(51) 



r-1 ' w l-t 

Notice that r c is only defined in the one-gap phase where t > t c , or equivalently 
S(t) > x- 111 terms of the new variables, we calculate the potential Uq at the circle^] 



N c " 2 ° g 5(t) - 1 



2r 



1, X 
-log 



1 



5(t) + 1 



(52) 



3 The potential is constant, otherwise the charge distribution would not be static. 
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and the interaction energy between the external particle and the charge distribution 



r(l — r) 



log 



5(t) + 1 
X + l 



+ Xlog 



4x 

(x + i)(*(t) + i) 



(53) 



Notice, that after the substitution Nf — > — £ these results are formally identical to 



the one-gap results of 27 where the bulk external particle had an opposite charge. 



However, there is an important difference in the definition of 5(t): in 27 S went to 



zero as the particle was approaching the unit circle, whereas in the present case S 
goes to infinity in this limit. 
The total energy becomes 



(2r 



N 2 

c 



r 
~2 



4 

log(l 



, X + l i, x- 



t 2 ) - 



(2r 



log 



1 r 2 A X 

log: 

-12 g (l + 5(t)Y 

X + l 1 , X - 1 r 2 

l l ° g W^-l + "2 hg W) 



X 



4 ~° 5(t) + 1 

giving the leading 0(N 2 ) term of loggjv /: c/(JVc)(*) i n ^ ne gapped phase: 



(54) 



(55) 



9N f ,U(N c )(t) 



X + l 
<J(t) + 1 

5(t) 

X 



(2N f -N c ) 2 /2 



x-i 

5(t) - 1 



JV c 2 /2 



;i-^) 



t r < t < 1 



(56) 



Notice that the results M9I) and (56) join smoothly for t — > t c , as 5(t) — >■ x- I n 



36 we proved that 0(iV c ) corrections to log g vanish, but there might be nontrivial 
contributions at (D(N®). Recall also that the result can be extended for unequal t 
and t by using the identity gN f ,u{N c )(t,t) = gN f ,u(N c )(^tt), i-e., by replacing t — > Vtt 
in the formulas above. 
As t — > 1, we get 



9N f ,U{N c ){t) 



x+i_y Nf - 



N c ) 2 /2 



X 



N 2 J2 



2N f N c -N 2 



(57) 



We observe that the exponent in the divergent factor (in the t — )■ 1 limit) is different 
from that of ( |49j ), A'" 2 has changed to 2NfN c — N 2 . This is good, because the power 
in the divergent factor then manifestly agrees with the dimension of the moduli space 
of U(N C ) SQCD, as expected (up to 0((N C ) ) corrections which vanish in this case). 

It is interesting to mention that the phase transition for gap formation is of third 
order, similarly as in 32 . This can be easily checked by calculating the derivatives of 
^°S9N f ,u(N c )(t) in (48) and (55) with respect to t and noting that the third derivative 



is discontinuous at t r 



3.2 Applying the log-gas approach in the SU(N C ) case 

We argue that the results of the previous sections hold for the SU (N c ) gauge group as 
well, if t = t, and up to subleading corrections in the Veneziano limit. This is expected, 
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because the definitions of the Hilbert series differ only by the extra condition for the 
total phase ^ a r a = 27m, which appears in the SU(N C ) case. When N c is large, 
small modifications of the configuration with 5r a ~ 1/N C are enough to make the 
total phase arbitrary while changing the contribution to the partition function by a 
sub leading amount]^] 

An analytic argument can be constructed as follows. We first separate the depen- 
dence of the total phase of the U(N C ) partition function: 

r d9 1 2j f 2w dr a ^ e n \ . . 

g NftUm (t,t) = / -x— H J" 2*5 YV.-0-27T* (58) 



2tt NJ 

N, 



xiA(,)i 2 n 



The overall phase 9 in the delta function can be eliminated by a small overall shift 
T a T a + 9/N c of the configuration. By rotational symmetry on the z-plane, this can 
be transferred to small displacements of the "external particles" at z = t and z = t, 
which leads to the relation 

g Nf , um (t,t) =f~ 9N f ,su(N c) (te- w / N %te^) (59) 

between the partition functions of the two gauge groups. 

Notice, however, that in the log-gas calculation above we set t = t. While this 
could be done for U(N C ) without loss of generality, it is a true extra constraint in the 
SU(N C ) case. For real t = t the effective Hamiltonian, given by the logarithm of the 



integrand on the right hand side of Eq. (59), becomes 



Nc AT , Q N c 

log \z a - z b \ + N f J2 log W -t\ + z2 1 _!| 2 + • '(S 60 ) 

l<a<b<N c a=l c a=l 1 a 1 

Taking into account the number of terms in the sums, the ^-independent term is 
O(N^) which makes it the leading one in the Veneziano limit, the linear term in 9 is 
0(N C ), and we dropped 0(N®) contributions. 

We shall now argue that the next-to-leading 0(N C ) correction vanishes in this case. 
This seems plausible, since due to symmetry with respect to the reflections along the 
real axis, gN f ,su(N c )(te l9 ^ Nc ) must be an even function of 9. More precise argument 
is as follows. In [36 ] we showed carefully that (minus two times) the minimal energy 
of the saddle point configurations gives the logarithm of the partition function up to 
0(N®) terms. Therefore, it is enough to show that the first variation of the saddle 



point energy due to the next-to- leading term oc 9 in Eq. (60) vanishes 



9 Actually, as we have verified numerically, this naive argument fails for t ^ t when t or t is close 
to one. Notice that we discuss here only the case where t — t where the log-gas approach is directly 
applicable. Then, as we shall see below, the U(N C ) and SU(N C ) Hilbert series indeed match in the 
Veneziano limit. 
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First, we note that the variation of the sum of energies of the leading O(N^) 
terms, which follows from small changes of r a in the saddle point configuration due 
to the additional term, is zero. This is the case because the saddle point energy 
is unchanged under any variation of r a by definition. Second, the leading saddle 
point configuration must be symmetric with respect to the real z-axis. Therefore, 
the possible 0(N C ) term obtained by evaluating on the saddle point configuration the 



next-to-leading term in the Hamiltonian (60), which is odd under the reflection with 
respect to real axis, vanishes. 



We conclude that for t = t the right hand side of Eq. ( 59 ) is independent of 9 up 
to 0(N®) corrections. Therefore 

log g N f,u(N e ){t) = logg NfjSU (N c )(t) + 0(N Q C ) (61) 

in the Veneziano limit. Therefore, the log-gas result holds up to (D(N®) corrections 
for both U(N C ) and SU (N c ) p*| However, we shall see numerically below that the con- 
vergence towards the leading term is considerably faster for U(N C ) than for SU(N C ). 



3.3 The GCBO approach 



Let us now compare our results with those in [16], which studied the behavior of 
the unrefined Hilbert series in the regime r = Nf/N c > 1 by using the Geronimo- 
Case-Borodin-Okounkov (GCBO) formula for Toeplitz determinants [37| . Their result 
reads 



9N f ,u(N c )(t, t) 

9N f ,SU(N c )(t,t) 



N c +1 



1 



[l-ti)r 



2 N? 



tly 2N c 

'rN c 

Nr. 



(62) 



F(N c ,r, tt)(t Nc +t 



N c + 1 



T(N c ,r,tt) 2 (tt) 



N c +1 



(63) 



where the numerators come from an expansion of the important determinant factor 
in the GCBO formula, and 



F{N c ,r,tt) = 2 F 1 {-{r-l)N c + l,rN c ;N c + 2;tt) 
F{N C , r, tt) = 2 F 1 {-{r - 1)N C + 1, r N c ; N c + 1; tt) . 



(64) 
(65) 



These expressions are expected to be good approximations for large N c and whenever 
tt(r — l) 2 is small. If r is fixed and tt is small, which is the limit we will be using 
here, N c does not need to be large. 

In the limit of large iV c , we expect the GCBO determinant factor, which gave the 
numerator in Eqs. (62), to contribute to the leading divergent behavior as t,t — > 1, 
so that the order of the divergence will correctly match with the dimension of the 



10 Actually, according to the GCBO result below the corrections in the gapless phase are exponen- 
tially suppressed. 
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moduli space. This cannot be however achieved with the above approximation which 
is restricted for small tt. 
Ref. 
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pointed out the need for good asymptotic formulas for the hypergeo- 
metric function. We study this by applying the saddle point method to the integral 
representation of the hypergeometric function 2 F\ and draw some lessons from the 
analysis. This analysis will show how the result (62) behaves in the Veneziano limit 



and for small tt, where it is expected to be applicable. 

The saddle point analysis is detailed in Appendix [DJ We find the following leading 



order results for the correction terms in Eqs. (62) and (63): 



rN c 

Nr. 



F(N c ,r,tt)(tt) 



N c /2 



rNr 



.N r . 



F(N c ,r, tt)(tt) 



z N c [F{s-)+^\ogtt\ 

cos(N c B + 4>) 



N c /2 

tt < (tt) 
tt > (tt) 



(66) 



c 

C ) 



where 



(tt) c 

F(s) 
s± 

s± 
B 



l/(2r- l) 2 

[r log s — (r — 1) log(s — 1) + (r — 1) log(l — sz)\ 



1 - tt ± ^(1 - tt) [l - (2r - l)Ht 
2rtt 



1 + 
1 
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1 -tt±iJ(l-tt) [(2r - l) 2 ti- 1 
2rtt 

(F(s + )-F(s-)) = ImF(s + ) 



tt < (tt) c 

tt > (tt) ( 



(67) 
(68) 

(69) 

(70) 
(71) 



and the order one phase is not determined at leading order (see Appendix D for 



the result at next-to- leading order). We see that there is a transition in the behavior 
of the correction terms, which takes place at tt = l/(2r — l) 2 . As the reader may 
have already noticed, the critical value matches with that of the log-gas analysis (see 
Eq. (|47|)), where a gap appeared in the charge distributions. In the saddle point 
analysis, the transition arises as the saddle points in the integral representation of 
the hypergeometric functions move from the real line into the complex plane, as can 
be seen from Eqs. (69) and (70), which give the saddle points in the two regions (see 



Appendix D for details) 



The saddle point result (66) for the correction terms in the GCBO approach 
behaves as follows. For small tt the factor -F(s_) + ^logtt is negative so that the 
corrections are highly suppressed. As ti approaches the critical value, F(s_) + | logtt 
approaches zero, and hence the corrections become large. For tt larger than the 
critical value, the corrections are 0(1) and oscillate wildly. This signals the fact that 



using only the first terms of the GCBO formula, which gave the results (62) and (63), 
is not sufficient to accurately describe the Hilbert series in this region. 

Finally, we compare the various results in Fig. [l] for t = t and for two choices of 
Nf and iV c . We also included the exact result, which can be quite easily evaluated 
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Figure 1: Comparison of the various approximations and exact results for (1 — 
t 2 ) Nf 9N f ,u(N c ){t) and (1 — t 2 ) N f gN f ,su(N c )(t) as a function of t. Top row: N c = 8 
and Nf = 16. Bottom row: N c = 8 and Nf = 32. Left column: linear scale. Right 
column: logarithmic scale. Thick (thin) lines are the results for the U(N C ) (SU(N C )) 
gauge group. Solid black lines are the exact results. Dashed blue lines give the GCBO 
result. The (single) dotdashed red line is the result in the Veneziano limit (given by 
the log-gas approach). 

numerically by using Toeplitz determinants of Appendix [A] up to N c ~ 10. The 
results are in line with our expectations. Generally, the GCBO result is a very good 
approximation at small t, but fails to reproduce the large t behavior. The leading 
term in the Veneziano limit, evaluated by mapping to log- gas, does not work as well 
for low t, but continues to be a reasonable approximation to the exact result all the 
way to t = 1. The GCBO approach seems to work better for SU(N C ) than for U(N C ), 
while for the log-gas result the opposite is true. 
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A Hilbert series as Toeplitz determinants 

In this Appendix we discuss the relation between the Hilbert series and Toeplitz 
determinants, pointed out in 16 . In the definition of the (refined) Hilbert series for 
the U(N C ) gauge group, 

1 Nc C 27T dr Nf Nc 1 

N c ] - ^ Jo 2tt (1 - tx t z a x ) (1 - ty iZa ) 

we may expand 

Nf oo 
f = \ (1 - ta^" 1 ) (1 - tyiZa) 

(recall that £ a = e JTa ). Then an easy computation (see, e.g., fl6]) shows that 

9N f ,u{N c )(t, i, x, y) = det T[f] = det (fi-j) ■ (74) 

V / t,j=l,...,N c 

While the Fourier coefficients /„ are quite complicated for the refined Hilbert series, 
in the unrefined case (setting all x, and y.j to unity) we find an explicit formula: 

f n (t,t) = (-ir(~ N /) 2 F 1 (N f + \n\,N f ,\n\ + l;tt)x { £ . (75) 



n 



Using this in Eq. (74) gives explicit results for the unrefined Hilbert series for small 
N c (and arbitrary Nf) which we shall use to numerically evaluate the Hilbert series 
exactly in some examples. 

Recall from Eq. @ that the Hilbert series for SU(N C ) could be written as a sum 
over integrals Ik otherwise similar to (72), but containing additional phase factors 
z\ • ■ ■ z% c . These modify the result (74) to read 



h(t, t, x, y) = det , (76) 

/ i,j=X,...,N c 



i.e., there is a shift by k in the indices of the Fourier coefficients. Again, the result (75) 
may be used to evaluate explicitly the integrals Ik(t,t), appearing in the expansion 
of the unrefined Hilbert series, for small iV c . 

Finally, we notice that the integrals for the unrefined series satisfy 

I k (t,t)=t k I k (tt) , (77) 



as can be verified by studying the determinant (|76j), or by studying the transformation 
integrals (|9| under rotations around the origin. In particular, the unrefined U(N C ) 
Hilbert series (which equals J ) depends on t and t only through the combination tt. 
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B Results for refined Hilbert series 



B.l U(N C ) with N f = N c + 2 



Here we prove Eq. (20 ). Starting from the general result ( 15 ), we can proceed similarly 



as for N f = N c + 1 (see Eq. (|T7j)): 

9N c+ 2,u(N c )(t,t,x,y) = ^(tt) |A| s A (x)s A (y) - (tif l sx(i)s x (y) (78) 

A X:£(X)=N c +l 

- (tt) lM sx(x)s x (y) . 

X:£(X)=N c +2 

The first term can be immediately summed, and the last term can be treated as in 
the Nf = N c + 1 calculation, while the middle term is more difficult to evaluate. 
It can however be summed by using a cofactor expansion. First, we notice that 
as £(X) — N c + 1, we can define a partition k such that k + l(iV c + 1) = A, i.e., 
Ki — Aj — 1 for all i — 1, . . . , -/V c + 1, and kn c +2 = 0. Now, as k runs over all partitions 
with £(k) < N c + 1, A covers all partitions with £(X) = N c + 1. By using the definition 
of s x , 

(79) 



A(x)sx(x) = det \ Xa 

/ i,j=l,...,N f 

Now the determinant may be processed by applying the cofactor expansion on its 
first row. Then the sub determinant of the kth term in the expansion is 



i,j=l,...,N c +l 



'N c +1 

n 



~(k) 



det[(x^ K(Nc+1) - i+1+1 " 



(80) 

i,j=l,...,N c +l 



-A(x^)s K (x^) 



X, 



where x^ k > is the sequence (xi, X2, ■ ■ ■ , x^-i, Xk+i, ■ ■ ■ , xjy f ), i-e, x with the element Xk 
removed. We find that 



Nr 



sx{x) 



fc=i 



1 A(x^) 
x\ A{x) 



s«(z<*>) ; (£(A) = iV c + l) . 



Inserting this in the second term of (78) we can sum over the partitions: 



£ (tt)^sx(x)s x (y) 

X:£(X)=N C +1 



(82) 



kl=l 



A(x) yj A(y) ^ 1 - ttx 
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where the sequence y^' is defined analogously to x^ above. Notice also that we were 
able to drop the constraint £(k) < N c + 1 since the modified sequences contain only 
N c + 1 variables. Combining with the other terms, and after some reorganization, we 
get 



9N C +2,U(N, 



N f 

)(t,t,x,y) = j [ 



X 



N f 



i _ at) N f - (tt) N < +i y — 



nSi(i - ustkik) n£i(i - ttxjm) 



fc,£= 



B.2 5t/(iV c ) with N f = N c + 1 



Here we prove Eq. ( 30 ) . The proof is quite similar to the one above for Nf = N c + 2 
with the U(N C ) gauge group. We start from Eq. (31) taking k > 0: 



I k = t kN ° t^s K+k[Nc) (x)^s K (y) . 

k:£(k)<N c 



Proceeding as above (see Eq. (80)) 



Nr 



m=l 



X k m +l A(X) K 



J5) 



Hence we get 



AT, 



EM) 



m,n=l 
AT, 



^ £ (-1) 



m+n 



m+n 



Y,(tt) ]K] s«(x im) )s«(y (n) ) (86) 



m,n=l 

A similar calculation for k < yields 



1 A(^ m )) 1 A(j/W) 
A(x) f n A(y) 

1 A(£( m )) 1 A(y^) 
xY 1 A(x) f n A{y) 11 i _ tfeMyW 



n 



(A; > 0) . 



j =f m c v r 1 A(£M) 1 A( ^ (n)) ft 

k 2^^ L > z Afft ~|fe|+l 11 



m,n=l 



A(x) A(y) Ali-ffe 



* y i 



(A; < 0) . 
(87) 
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Inserting the result in the definition (24) the sum over k can be done, giving 

9N c +i,su(N c )(t,t,x,y) ( 

y r })m+n ( 1 i 1 ,U AgM) 1 A(yW) 

" ' ll-^™ 7^ A(i) #1+1 A(y) 



i,j=i 1 '< yj 



= TT — ^- 

1,3=1 iy 3 

x ^ (i - (gyyg^)) n£(i - g*go ngia - 

j^i^C 1 - t N */x k )(\ - t Nc /ve)(l - ttXkVe) Ui^i - £ fc) Uj^eiVj ~ Vt) ' 

C Asymptotics of the unrefined Hilbert series 

In this Appendix we analyze the asymptotics of the unrefined Hilbert series and the 
singularity structures of the consequent sums. First, we recall that 

oo 

9N f ,u(Nc){t,t) = 5>t) s G(s,0) (89) 

s=0 

oo oo —1 oo 

9N f ,su(N c )(t,t) = ]T]Tt feM (ttrG( S ,A;)+ E Y,W M (i£) a G(s,\k\) , (90) 



k=0 s=0 fe=-oo s=0 



where 



G(s,k) = ^2 d x+k (N c )d\ (91) 

A:|A|=s, £(\)<M 

M = min (N c ,N f ) . (92) 

Therefore, we need to find the asymptotics of G(s, k) for large s and k. We do this 
by a "brute force" calculation using the standard formula 

- LJ - Mm) 

me A v ' 

where the box m of the Young diagram A lies on the z(m)th row and the j(m)th 
column, and h(m) is the hook length associated to m: 

h(m) = \ j(m) - i{m) + A- (m) - j(m) + 1 (94) 

with A' being the conjugate partition, i.e., \[ is the length of the ith row of A. 
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Let us abbreviate Nf = N. We shall divide the diagram A into M rectangular 
blocks, denoted by A^, and defined by 

me \W^K {m) = K , K=1,...,M, (95) 

i.e., the Kth block consists of rows that have the length K. Then we have 

M 



K=l 



dx = Y[d x(K) (96) 

'(m) — 
h(m) 



dxw = 11 • (97) 



m£X( K ) 

Since the blocks are rectangular, it is straightforward to write d X (K) in terms of Euler 
gamma functions. At this point, let us also reintroduce the variable /cj^j] We get 

X K +k K N ^ _ . 

<wo)<*> = n u Xj+ k-ilK-j+i m 

i=X K+1 +k+l j=l ■> 
A 



T(N + x K + k- 3 + i) r(A,- - Aa: + a: - j + 1) 
r(A + x K+1 + k-j + i) r{Xj - \ K+1 + K-j + i) 



n 



for 1< A < M and 



A? 



_ TJ T(N + X M + k - j + 1) r(Aj - A M + M - j + 1) 
<W(M))(M) - 11 T{N _ j + 1) T{Xj + k + M-j + l) • ( " } 

In order to find G(s,k) at large s = |A| = J^A, an d & we need to analyze these 
expressions in the limit where Aj — > oo with their ratios fixed, and also k — > oo. We 
shall not take s/k fixed, but rather consider limits where either of them is taken to 
infinity independently. The analysis can be done by using the asymptotics of the 
gamma functions, i.e., the Stirling formula 

logT(A + x) = X log A - A + (x - log A + l - log(27r) + O Q J , (100) 

where we neglect the constant term, since we are not trying to calculate the propor- 
tionality constant of the final scaling results. We also notice immediately that the 

n We shall consider the partition A + k(M), i.e., not A + k(N c ). Whenever N c ^ M {i.e., M = 
Nf < N c ) only the k = term contributes in the sums above, and Sx-i-k(N e ) (^0 = 0- 
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term ~ —A does not contribute due to cancellations. Denoting f(x) = xlogx, we get 



d 



(101) 



(X+k(M))( K ) 

r K K K-l K 

~ exp J2f(^K+k)-J2f( X K + i+k) + J]/(A i -A K )-J]/(A J -A K+1 ) 

L j=i j=i j=i j=\ 

r- K K 

x exp ^(N - j + 1/2) log(A^ + A;) - £(JV - j + 1/2) log(A x+1 + A;) 
L j= i j=1 

r-K-l K 

x exp ^(K-j + 1/2) log(A J -A i ,) - ^(K - j + 1/2) log{Xj-X K+1 ) 

for 1 < < M. We left out all factors that are finite as s — > oo (with all Aj/Aj fixed) 
or — >■ oo. Notice also that the sum of the log-terms (last two rows in ( 101 )) is finite 
in both limits if the j — K term in the last sum is removed. Therefore, we shall only 
keep this divergent term below. For K = M we find 



d 



exp 



M-l 



M 



(\+k(M)Y M ) 
M 

3=1 3=1 3=1 

M 

x exp 

L 3=1 

rM-1 

x exp 

L 3=1 3=1 



(102) 



+ 1/2) log(A M + fc) 

3=1 
M-l 

^ (M - j + 1/2) log(A 3 - -A M ) - £(M - j + 1/2) log^+fc) 



Let us then insert the results in the formula (96). It is not difficult to check that all 
terms involving / cancel. We are left with 



d 



X+k(M) 



exp 



M-l 



M 



M-l 



l - log{X K -X K+1 ) + jf^JV ~ J + 1/2) log(A M +A:) (103) 

K=l j=l 

M 

+ (M - j + 1/2) log(A, -A M ) - £(M - j + 1/2) log(X 3 + k) 

3=1 3=1 

To continue, we write an expression which correctly interpolates between the leading 
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asymptotics as either s = |A| — )■ oo or k — > oo first r^j We get formally 
d 



exp 



M-l 



M 



log s + j^(iV - j + 1/2) log(s+Cfc) (104) 

*sr=i i=i 

M-l M 

+ J2(M-J + 1/2) log s - £(M - j + 1/2) log(s + Ck) 

3=1 j=i 
„ ( s+C7A .)M(iV-M) s M(M-l)/2 ^ 

where C can be treated as a constant. The remaining task is to sum over the partitions 
in the definition (91 ). As the scaling result (104) is independent of Aj, this amounts to 
calculating the number of partitions that enter the sum. We sum over Ai > A2 > • • • > 
Am with one constraint s = ^ Aj, so the number of terms is ~ s x . Therefore 



13 



G{s,k) 



~ s M l d\+k(M\d\ 



*A+fc(M)"A+0(M) 



(105) 



\M(N—M) M(M-l)/2 M(N-M) M(M-l)/2 



s"' ~(s+Ck) ^' "*'s-" v - ••••-«-••*• 

( s + CA .)M(iV-M) s MiV-l 



The main results are given in the text where we assume that Nf > N c so that 
M = N c (and iV = Nf). However, the analysis also works when Nf < N c and 
M = Nf, even though only the k = case appears in the sums (90) in this region. 
Indeed, we find 

G(s,0) ~ s^f- 1 (106) 



and 



9N f ,U(N c )(t,t) = 9N f ,SU(N c )(t,t) ~ ^2(tt) s S N f 1 



1 



s=0 



(107) 



for Nf < N c as tt — > 1, which is the expected result (compare with (16) with 
y 3 = 1 for all i,j). 



D Saddle point analysis of the GCBO result 



In this Appendix we shall calculate the leading asymptotics of the results (62 ) and (63 ) 
in the Veneziano limit. We start from the standard integral formula 

2 F 1 (a, b; c; z) = ^ [ dss b ~\l - 8 ) M {1 - zs)~ a . (108) 

r(6jr(c- b) j 

12 A more precise way would be to simply treat the limits separately. Then we would not need to 
introduce the 0(1) functions C which depend in a complicated way on Aj/s and will cancel in all 
results. 

13 The result ([105} covers all nontrivial cases, i.e., it holds for k = or Nf > N c . For k ^ and 
Nf < N c the result is not applicable since we used the formula for d\ + k(M) instead of d\ + k(N c )i but 
in this case actually G(s, k) = 0. 
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However, in our case b > c so the integral will not converge at t = 1. We then regulate 
the t = 1 singularity in a standard manner: we deform the integration contour and 
write 



2 Fi(a,b; c; z) 



r( c ) 



1 



r(6)r(c - b) 1 - e 2 ™( c - fc ) 



6-1/ 



zs) 



(109) 



where the integration contour starts from the origin, circles around s — 1 in the coun- 
terclockwise direction, and returns to s = via a different branch of the integrand. 
We can contract C to two pieces of a real line and a circle around s = 1 having radius 
e. In the domain Re(c) > Re(6) the contribution on the circle vanishes as e — > and 



we recover the original representation (108), which proves the formula (109). 





Figure 2: An illustration for the saddle point method on the complex s-plane. Left: 
saddle points at real values of s. Right: saddle points at complex s. The contours in 
the background are at constant values of Re F(s) with lighter shades denoting smaller 
values. The filled circles are singularities of F(s), while the open circles are the saddle 
points. The integration contour C is the thick curve that starts at s = 0, goes around 
the singularity at s = 1 through one (left) or both (right) of the saddle points, and 
returns to s = 0. 

Now the integral representation converges in the domain of interest to us and we 
can proceed with the saddle point analysis. Let us first rewrite the representation as 



2 F 1 (a,b; c; z) 
I 



r( c )r(i + 6-c) 



i 



r(6) 



(110) 
(111) 



We will restrict to the leading result (next-to-leading corrections will be considered 
below) so we take a = — (r — 1)N C , b — 1 = rN c and c — b — 1 = — (r — l)iV c which 
covers both functions in (64). Also recall that z = tt lies between zero and one. The 



crucial function is the logarithm of the integrand: 

F(s) = [r log s — (r — 1) log(s — 1) + (r — 1) log(l — sz)] 



112) 
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where we left out the common factor N c . Saddle points are found at the zeroes of 

r r — 1 z(r — 1) — rzs 2 + s(l + 2rz — z) — r 
F ( s > = - S -J^- = s(s-l)(l-sz) (113) 



and are given by 



, ± = 1+ i-*±^(i--)[i-(fr-i)'--]. (114) 

2rz 



We observe that there are two possibilities, separated by the critical value z = tt = 
l/(2r — l) 2 . We immediately notice that this matches with the critical value of t 
where the gap appears in the Coulomb gas approach. Let us analyze the two cases 
separately. 

• For z = tt < l/(2r — l) 2 there are two real roots. It is easy to see that they lie 
between s — 1 and s = 1/z. Since F'(s) — > — oo when s approaches one from 
above, F"(s-) > and F"(s+) < 0. We choose C such that it passes through 
s = S- perpendicularly to the real axis. The leading behavior is then 

/ ~ e "oF(,_) ) (115) 

where the proportionality constant is real and positive. We did not find a simple 
formula for F(s_). Using this result we find that 

rN c \ , at /i ( rN, 



J(N c ,r,z)z N ^ 2 ~ ( A \ )j(N c ,r,z)z N ^ 2 ~ e ^[F(--)+|iog*] ^ (n6) 



N C J v ' ' \N C + 1 

where the binomial coefficients and gamma functions cancel at leading order. 
It is straightforward to show that the combination F(s-) + (log2;)/2 is negative 
for z < l/(2r — l) 2 (and approaches zero as z — » l/(2r — l) 2 ). Therefore, the 
corrections are exponentially suppressed in N c in the U(N C ) case. For SU(N C ) 
gauge symmetry, corrections are suppressed when t = t. For t ^ t there are 
exponentially growing corrections for z sufficiently close to l/(2r — l) 2 , as the 
result for J 7 in Eq. (116) will be multiplied by 

t \N c /2 /f\ N °/ 2 



ti +u> • (117) 

For z = tt > l/(2r — l) 2 there are two complex roots: 

l-z± iJ(l - z) \(2r - l) 2 z - 11 . . 

s± = 1 + ^ ^ ' . 118 

2rz 

Let us denote F(s±) — A ± iB. A straightforward calculation yields A = 
— (logz)/2. We find that all exponential factors cancel: 



rN, 



^(N c ,r,z)z N ^ 2 ^ ^ N ^T(N c ,r,z)z N ^ 2 ^ CO s(N c B + <P) . (119) 

The vividly oscillating 0(1) result suggests that the approximation of Eqs. (62) 
and (63) is breaking down. 
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Finally, let us consider higher order corrections to the saddle point results. The 
exact form of the logarithm of the integrand in the integral representation of the 
hypergeometric function is 



F(s) = N c F(s) + G(s) 



(120) 



where G(s) differs for the two expressions of Eq. (64): the parameters a = — (r — 



1)N C + 1 and b = rN c of the hypergeometric functions are the same in both cases, 
but c is different. Therefore, we define 



G(s) 



log(s - 1) - log(s) - log(l - zs) = Gi(s) 
- log(s) - log(l - zs) = G 2 (s) 



c = N c + 2 
c = N c + l 



'121) 



We shall denote the "corrected" zeroes of F'(s) as s±. The full Gaussian contribution 
to the saddle point result, 



1 



2ttF"(S_ 
1 

2nF"(s_ 



for real roots 
1 



2nF" '{S. 



F(S_) 



(122) 
(123) 



2Re 



27rF"(s + ) 



3 F(S+) 



for complex roots 



is accurate up to 0(1/N C ) corrections. Expanding in 1/N C we find, after some calcu- 
lations, 



1 



t/2ttN c F»(s-) 
1 



N c F{s-)+G(s-) 



l + O 



2Re 



y/2nN c F"(s+) 
The correction terms become 



1 



l + O 



tt<l/(2r-l) 2 (124) 



tt > l/(2r- 1) 



rN c 
N c + 1 



rN c 



T(N c ,r, z)z 



iV c /2 



T(N c ,r,z)z N ^ 2 



^27cN c F"(s.)r-l 



(125) 



x exp 



iV c F( S _) + ^logz + G 1 ( S _) 



^2TiN c F"(s-)r-l 



(126) 



x exp 



N 

N c F(s^) + ^ log z + G 2 (sS / 
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in the region of real saddle points, and 

U+iK«'^ - 7^Sm\^' MN ' B+M (127) 

in the region of complex saddle points, up to 1/N C corrections. Here B = InxF(s + ) 
as above, 

0, = -iaxgF"(s+) + ImGi(s + ) , (128) 

and we used the results 

ReG 1 (5 + )=0; Re G 2 (s+) = - log — -. (129) 

2 (r — — z) 
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